Macromolecular unfolding properties in presence of compatible solutes 
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We present Molecular Dynamics simulations of Chymotrypsin inhibitor II and PEO in presence of 
compatible solutes. Our results indicate that the native compact state of the studied macromolecules 
is stabilized in presence of hydroxyectoine. We are able to explain the corresponding mechanism 
by a variation of the solvent properties for higher hydroxyectoine concentrations. Our results are 
validated by detailed free energy calculations. 
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1. INTRODUCTION 



The cell metabolism of microbacteria at extreme tem- 
peratures and dry environments is protected by the pres- 
ence of compatible solutes which are also called os- 
molytes, co-solvents or extremolytes [iHll- It has been 
found that compatible solutes are strongly hygroscopic 
and prevent denaturation processes of the cell proteins by 
increasing the melting temperature [H 0, @] ■ One impor- 
tant osmolyte which is often used in industrial products 
is hydroxyectoine 0, Especially in the last years, sev- 
eral studies have been published to investigate the gen- 
eral function of osmolytes. Most of the theories which 
aim to describe the stabilization mechanism accurately 
can be divided into direct and indirect interactions 9- 
|. In the context of extremolyte protein stabilization, 
we will briefly discuss in the following the most common 
ones. 

The framework of the preferential exclusion model as an 
indirect mechanism is established by the consideration 
of a protein dissolved in aqueous solution that contains 
compatible solutes Within this model it is as- 

sumed that the chemical potential of the water at the 
proteins surface is perturbed by the presence of the co- 
solute. The corresponding thermodynamic interactions 
[l3| lead to deviations from the bulk molar concentra- 
tion of the co-solvent compared to the local environment 
around the protein. These deviations can be related to 
a preferential binding or a preferential exclusion of the 
co-solute. 

The dynamic exchange of extremolytes in the local sur- 
rounding of the proteins surface result in free energy dif- 
ferences due to exclusion. These differences are compen- 
sated by a significant increase of the water molecules in 
the local environment of the protein. Preferential hydra- 
tion and preferential exclusion are therefore synergetic ef- 
fects which appear simultaneously. According to the orig- 
inal idea, the stabilizing effect on the protein is mainly 
caused by the excess solvent molecules whose interac- 
tions restrict the protein to its native state by increas- 



structure in presence of ectoine, trimethylamine-N-oxide, 
amino acids and sugars has been also observed in exper- 
iments Isrl - liol l . In addition, the preferential exclusion 
model has been used to understand the effects of mono- 
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ing the melting temperature [15l . Il6l | . A varying melting 



temperature as well as a strong influence on the water 



layer structure formation 14 1. 
Molecular Dynamics simulations have indicated 
that the preferential exclusion model is relevant for a 
specific class of proteins. It has been found that the co- 
solvents do not directly interact with the proteins surface 
in agreement to the theory, but slow down the diffusion of 
solvent molecules in the local environment of the protein 
[l5| . It has been concluded that this restriction is re- 
sponsible for the change of the melting temperature and 
therefore the sustainment of the kinetic stability of the 
protein. This has been in particular investigated for the 
Chymotrypsin inhibitor II (CI2) ^16']. Other proteins like 
Met-Enkephalin show deviations to this model which re- 
sult in the authors conclusion, that the preferential exclu- 
sion model is only relevant for proteins with well-defined 
and pronounced hydration shells 16|. Macromolecules 
which are less well-hydrated, e. g. Met-Enkephalin are 
only slightly influenced in their unfolding processes such 
that the melting temperature remains unchanged. 
A different mechanism in contrast to the preferential 
exclusion assumption is proposed by the transfer free 
energy model [l3, [l^l which relates the increased sta- 
bility of the protein to a solvent quality decrease in 
presence of compatible solutes. This leads to favorable 
intra-molecular hydrogen bonds compared to unfavorable 
solvent-protein interactions. The excess of intramolecu- 
lar hydrogen bonds results in an increased stability of 
the native structure which is validated by a shifting of 
the equilibrium constant to the folded state. 
In this paper we focus on the properties of CI2 and PEO 
(polyethylene oxide) in presence of compatible solutes. 
Our results indicate that preferential hydration is absent 
for these molecules although an increased stability of the 
native state can be observed. We have found that the 
free energy landscape is significantly changed in presence 
of hydroxyectoine such that unfolding becomes energet- 
ically less favorable. In addition we present numerical 
findings which validate that changes of the solvent prop- 
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erties are mainly responsible for the observed stabiliza- 
tion behavior. 

The paper is organized as follows. In the next sections we 
introduce the theoretical background and the numerical 
details. The results are discussed in the fourth section. 
Our conclusions and a brief summary are presented in 
the fifth section. 



2. THEORETICAL BACKGROUND 



Structural properties of proteins an d polym ers can be 



described by the end-to-end radius 
to-end vector is given by 

re = rM - r 1 



41-41. The end- 



(1) 



where rjy/i denotes the position of the A^-th and the first 
monomer. The end-to-end radius obeys a scaling behav- 
ior which is expressed by 



2v 



(2) 



where v denotes the excluded- volume parameter [41|-|43j 
which is dependent on the solvent properties 4^ . Hence 
a good solvent for the macromolecule results in a swollen 
conformation whereas a bad solvent leads to a shrinkage. 
A general way to calculate free energy differences for 
macromolecular conformations or solvation properties in 
computer simulations is thermodynamic integration (45| . 
The differences can be computed by the coupling of the 
Hamiltonian H to a perturbation parameter A which cor- 
responds to a well-defined state [46| . This parameter reg- 
ulates the strength of the conservative interactions from 
fully developed to absence. The free energy difference 
between state A and state B is defined by 

/\Fab = F{A) - F{B) = ^^dX 



dX 



(3) 



Xa 



where A is in the parameter range from to 1. This 
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FIG. 1: Thermodynamic cycle for the calculation of solvation 
free energy differences. 



equation can now be used to calculate free solvation en- 
ergies 4^ 48 1 where A = represents the dummy state 
without intra- and intermolecular interactions. The cor- 
responding thermodynamic cycle is shown in Fig. [TJ The 
solvation free energy is then defined by 

AF = AF3 - AFi (4) 

where AF2 vanishes due to non-existing interactions. 

3. NUMERICAL DETAILS 

We have performed Molecular Dynamics simulations 
in explicit SPC/E water model solvent [H^l with the soft- 
ware package GROMACS [HlUHsj and variations of the 
the force field GROMOS96 54]. We created the chem- 
ical structure of hydroxyectoine ((4S,5S)-2-methyl-5- 
hydroxy-l,4,5,6-tetrahydropyrimidine-4-carboxylic acid) 
as it was in detail described in Ref. @. 
The Molecular Dynamics simulations of Chymotrypsin 
inhibitor II (CI2) have been carried out in a cubic simu- 
lation box with periodic boundary conditions. The box 
consisted of (5.76637 x 5.76637 x 5.76637) nm^ filled with 
5498 SPC/E water molecules. The presence of hydrox- 
yectoine was simulated by a 0.087 mol/L solution which 
corresponds to 10 molecules. CI2 was taken from the pdb 
entry lYPB 

The polyethylene oxide (PEO) chains consisted of 31 
monomers where the to polo gy of the chain was modeled 
by the PRODRG server [56!| . The box has a side length of 
(5.32547 X 5.32547 x 5.32547) nm^ filled with 4944 water 
molecules. The presence of hydroxyectoine was simulated 
for several concentrations between to 0.11 mol/L which 
corresponds to ten molecules. 

Electrostatic interactions have been calculated by the 
Particle Mesh Ewald sum 57| . The time step was St = 2 



fs and the temperature was kept constant by a Nose- 
Hoover thermostat 45] at 300 K for the PEO molecule 
and 400 K for the Chymotrypsin inhibitor II. All bonds 
have been constrained by the LINCS algorithm 58]. Af- 
ter minimizing the energy, we performed 1 ns of equili- 
bration followed by a 10 ns simulation sampling run. 
Solvation free energy calculations for ectoine ((S)-2- 
methyl-l,4,5,6-tetrahydropyrimidine-4-carboxylic acid) 
have been performed by independent simulation runs for 
parameter values of A 6 { 0, 0.05, 0.1, 0.15, 0.2, 0.25, 
0.3, 0.35, 0.4, 0.45, 0.5, 0.55, 0.6, 0.65, 0.7, 0.75, 0.8, 
0.85, 0.9, 0.95, 1}. Each configuration has been equi- 
librated for 250 ps and the computation of free energy 
values has been achieved in a 900 ps simulation run for 
each A. 

To avoid singularities, the ordinary potentials have been 
replaced by the following soft core interaction 



(l-A)V^'^ (a(T'^A + rS)s 
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(aCT6(l-A) + r^)« 



(5) 



where V{r) denotes the ordmary electrostatic or 
Lennard- Jones potential in the GROMOS force field in 
the A-state (A = 0) and in the B-state (A = 1) [59]. The 
parameter a has been chosen to 0.7 and the radius of 
interaction value a was 0.3. 

The free energy landscape for the PEO chain has been 
calculated by the metadynamics method presented in 
Ref. [60j . A history dependent biasing potential in form 
of Gaussian hills was applied to the molecule during the 
simulation which helps to overcome energetic barriers. 
Details of the method can be found elsewhere [60l - l62j |. 
The simulations have been conducted by the program 
plug- in PLUMED [il] . The Gaussian hills were set each 
2 ps with a height of 0.1 kJ/mol and a width of 0.25 nm. 



4. RESULTS AND DISCUSSION 

4.1 Chymotrypsin inhibitor II 

We start this section by discussing the results for Chy- 
motrypsin inhibitor II in presence of hydroxyectoine. A 
snapshot of a typical configuration for CI2 in a 0.087 mo- 
lar solution of hydroxyectoine is shown in Fig. [2] At a 
first glance it can be seen that direct interactions of the 
co-solvent with the protein are nearly absent. For a more 
quantitative description, we have calculated the average 
number of intermolecular hydrogen bonds between sol- 
vent and protein nj^'" and between the residues of the 
protein n^^^ for two different concentrations of hydrox- 
yectoine shown in Tab. HI 

It is clearly visible that the values for nu are nearly iden- 
tical for both concentrations. This becomes in particular 
true by regarding the values for the intermolecular hydro- 
gen bond density, which is defined by p^"'" =< n|^''"/E > 
with the total solvent accessible surface area E. The val- 
ues are p^°^'" — 1.93 nm^^ in presence and p^°^'" — 1.94 
nm~^ in absence of hydroxyectoine. 
Furthermore only slight variations of the intramolecular 
hydrogen bonds n^^^^ can be observed. It can be con- 
cluded that the results do not agree with the preferential 
exclusion [l3| due to identical values for the hydrogen 
bond densities as well as with the standard transfer free 
energy model (l^ due to the negligible differences in the 
number of intramolecular hydrogen bonds. 
However, the stabilization effect of hydroxyectoine is 
clearly evident by analyzing the average end-to-end radii 
Te whose values are also shown in Tab. HI It can be seen 
that the protein is in a more compact conformation in 
presence of hydroxyectoine which is in agreement to the 
results published in Refs. [H, Finally we have cal- 
culated the solvent accessible hydrophilic Sp''*' and hy- 
drophobic surface area YP^°^ , It comes out that the hy- 
drophilic surface area decreases while the hydrophobic 



TABLE I: Average number of intermolecular hydrogen bonds 
n|f'" between the solvent and CI2 and intramolecular hy- 
drogen bonds n%''^ between the residues of CI2, the average 
end-to-end radius Ve and the solvent accessible surface hy- 
drophilic/hydrophobic area 

^phii/phob f^j. concentrations 

c of hydroxyectoine. 



c [mol/L] 










^phob 


0.087 


117.84 ±0.01 


30.88 ±0.05 


1.45 ± 0.01 


37.7 


22.8 





115.58 ±0.05 


33.48 ± 0.04 


1.54 ±0.01 


38.1 


22.0 



area increases in a hydroxyectoine/water mixture. These 
at a first glance astonishing results can be directly re- 
lated to a change in the solvent properties and will be 
discussed in the last section. 

In addition we have studied the number of hydrogen 
bonds between CI2 and hydroxyectoine. We have found 
an average value of 1.09 ± 0.01 hydrogen bonds between 
both solutes. This result is in agreement to a recent 
publication, which has reported the occurrence of direct 
binding between hydroxyectoine and lipid headgroups for 
high molar concentrations 6J|. Nevertheless, regarding 
the large number of hydroxyectoine molecules, it can 
be concluded that preferential binding of hydroxyectoine 
which has been found as important |21| for zwitterionic 
co-solvents is a minor effect for protein stabilization. This 
becomes obvious by regarding the correlation coefficient 
between the end-to-end radius and the number of direct 
bonds with hydroxyectoine which is given by a negative 
value of —0.012. Hence it can be concluded that the sta- 
bilization of a compact state in presence of hydroxyec- 
toine is mainly caused by indirect solvent interactions. 



4.2 Polyethylene oxide 

To study the stabilization process in more detail, we 
have simulated a PEO chain in aqueous solution for var- 
ious hydroxyectoine concentrations. PEO consists of po- 
lar and apolar atoms {jyha lYphob _ qj) ^^^^ ^-^^^ ^j^g 

stabilizing effects of co-solvents can be investigated by a 
chemically simple PEO molecule instead of complex pro- 
teins. 

We start this investigation by calculating the end-to-end 
radius re for PEO in absence of hydroxyectoine and for 
a concentration of 0.11 mol/L. The average end-to-end 
radius is r,, — 1.17±0.01 nm in absence of hydroxyec- 
toine and Tf. ~ 0.91 ± 0.01 nm for a 0.11 molar con- 
centration. Thus PEO shows a significant shrinkage of 
about 22% for the end-to-end distance. To investigate 
the influence of direct interactions between PEO, hy- 
droxyectoine and water molecules, we have calculated the 
number of intermolecular hydrogen bonds with water. It 
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FIG. 2: Simulation snapshot of Chymotrypsin inhibitor II in 
presence of hydroxyectoine. 



is nearly constant for all hydroxyectoine concentrations 
with a value of 7^^°'"" — 5.38 ± 0.02. The same is true for 
hydrogen bonds between hydroxyectoine and PEO with 
^PEO ^ 0.07 ±0.02. Thus it is obvious that direct inter- 
actions between PEO and hydroxyectoine are absent. 
The influence of varying hydroxyectoine concentrations 
on the end-to-end radius can be observed in Fig. |3l It 
becomes obvious that similar to CI2 the PEO chain is 
more compact in presence of hydroxyectoine. This effect 
is drastically pronounced above a concentration of 0.09 
mol/L. 

To study the swelling process energetically, we have cal- 
culated the free energy landscape for the end-to-end ra- 
dius via metadynamics simulations. A nearly identical 
method has been applied in a recent publication [s^ . The 
landscape in presence (c=0.11 mol/L) and in absence of 
hydroxyectoine is presented in Fig. U An increase of the 
free energy can be observed in presence of compatible 
solutes for larger values of the end-to-end radius. For a 
pure solution, re values between 0.55 to 1.65 nm with an 
energy difference of AF « 0.3 kcal/mol are energetically 
comparable. In presence of hydroxyectoine, identical dif- 
ferences can be observed between ~ 0.55 and 0.85 nm. 
Values larger than 1 nm are energetically unfavorable in 
presence of the extremolyte. Combined with the above 
presented results, it can be concluded that the compact 
structure is stabilized due to an increase of the free en- 
ergy barriers in presence of hydroxyectoine. 



1.2 




0.9 I ' ' ' ' ' 1 

0.02 0.04 0.06 0.08 0.1 0.12 

c [mol/L] 

FIG. 3: End-to-end radius for the PEO chain in presence 
of varying hydroxyectoine concentrations c. 
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FIG. 4: Free energy landscape for the end-to-end radius of 
PEO for a 0.11 molar hydroxyectoine solution (triangles) and 
and in absence of hydroxyectoine (circles) . 



4.3 Solvent properties 

To study this effect in more detail, we have calculated 
the free solvation energy of ectoine, which is a zwitte- 
rionic molecule for various concentrations of hydroxyec- 
toine. The results for the free solvation energy differences 
to a pure solution AF° are shown in Fig. [5] It is worth to 
notice that the free solvation energy of ectoine is largely 
exothermic with a value of Af^ = —89. 97 ±0.1 kcal/mol. 
Thus it can be concluded that the solubility of ectoine is 
decreased for larger concentrations of the co-solvent. The 
linear dependence is obvious. Direct interactions between 
the co-solvent and ectoine are absent. Hence it becomes 
obvious that the main mechanism is indirectly mediated 
by the solvent such that the free solvation energy for po- 
lar residues or atoms is significantly increased in presence 
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FIG. 5: Free solvation energy differences of ectoine for var- 
ious concentrations of hydroxyectoine compared to the pure 
solution. 
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FIG. 6: Excess free energy difference to the pure solution AF'' 
for various concentrations of hydroxyectoine. Inset: Corre- 
sponding excess entropies AS" . 



of hydroxyectoine. 

To study the consequences for the solvent in detail, we 
have performed simulations with varying concentrations 
of hydroxyectoine in absence of macromolecules or other 
solutes. 

The excess free energy yields an estimate for the devia- 
tion of a real solution to the ideal state. It can be com- 
puted by the following equation 



F"^ = F-F,d^^kBT log <e 



-V/knT 



> 



(6) 



where V denotes the potential energy of the mixture and 
ksT the thermal energy [4^. The excess free energy 
difference of a hydroxyectoine/water mixture to a pure 
solution of water is given by Ai^° = F'^^ — F'^^''^ where 
pex,o jg ^Yie excess free energy of the pure water. The 
results are shown in Fig. [6] The increase of the excess 
free energy difference for larger hydroxyectoine concen- 
trations is remarkable. 

Furthermore the excess entropy difference is given by 
AS° = (AC/° - AF°)/T with AC/" = U"^ - C7^=^'° where 
U is the total energy of the system. In agreement to the 
excess free energy differences, decreased entropies can be 
observed for larger concentrations of hydroxyectoine. It 
can be concluded that due to the decreased entropy dif- 
ferences the solvent is locally more ordered in presence 
of hydroxyectoine in agreement to recent results 0] . 

The effects of hydroxyectoine on the solvent properties 
can be also observed by tlie_ investigation of the dipolar 
autocorrelation function 



< /2(t)/2(io) >~ exp(t/T) 
T averaged over all water 



with the correlation time 
molecules. The results are presented in Fig. [71 Although 
the effects are small, the linear dependence of the con- 
centration on the relaxation times as well as the aver- 
aging effect over all molecules emphasizes the relative 
influence of hydroxyectoine. It is obvious that the pres- 
ence of a high molar concentration of hydroxyectoine is 
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FIG. 7: Normalized dipole autocorrelation function for wa- 
ter molecules in varying hydroxyectoine concentrations. In- 
set: Relaxation times r for the different hydroxyectoine 
concentrations. The values can be fitted by regression to 
TO = 4.82 ± 0.01 ps with a slope of 1.77 ± 0.05 ps(mol/L)-\ 



directly related to increased relaxation times. Hence it 
can be concluded that the structural reorganization of 
water is diminished in presence of hydroxyectoine which 
validates the observed excess entropy values shown in 
Fig. [51 Similar effects can be also observed in the con- 
text of hydrophobic hydration [66| and the influence of 
disaccharides on water dynamics [67| . 
The decrease in entropy is also supported by the results 
shown in Fig. [8l where the average number of hydrogen 
bonds for a single water molecule is presented. Although 
the variations are small, it can be assumed that the gen- 
eralized net effect of a smaller number of hydrogen bonds 
for all water molecules compared to the zero hydroxyec- 
toine concentration is significant. The above discussed 
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FIG. 8: Average number of hydrogen bonds of a single wa- 
ter molecule with other water molecules for varying hydrox- 
yectoine concentrations. Inset: Corresponding lifetimes for 
hydrogen bonds. 
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FIG. 9: Diffusion coefficient of water molecules for different 
hydroxyectoine concentrations. 



increased relaxation times are additionally included in 
the hydrogen bond lifetimes which show an increase of 
0.03 ps. 

All the above presented results are related to a signifi- 
cant variation of the entropic contributions and the local 
ordering of the solution resulting in an increase of the 
free energy. 

A restriction in molecular motion becomes also obvi- 
ous in Fig. [9] where the average diffusion constant D for 
a water molecule in presence of varying hydroxyectoine 
concentrations is shown. It is evident that D is smaller in 
presence of larger hydroxyectoine concentrations. It has 
to be noticed that this result is in agreement to Ref. [3] , 
where the slowdown of water molecules for higher ectoine 
concentrations has been additionally observed. Thus we 
assume that the presence of hydroxyectoine results in 



a decreased water diffusivity in agreement to Ref. |l5| . 
Thus we conclude that hydroxyectoine changes the prop- 
erties of the solvent such that unfolding of proteins be- 
comes less favorable. 



5. SUMMARY AND CONCLUSION 

We have investigated the influence of compatible so- 
lutes on the structural properties of CI2 and PEO. It has 
been found that in presence of hydroxyectoine the protein 
as well as the polymer show a significant preservation of 
the compact structure. Energetic stabilization of com- 
pact states due to an increase of the free energy barriers 
has been observed by analyzing the results of metady- 
namics simulations. These findings are supported by an 
increase of the free solvation energy in presence of hy- 
droxyectoine for polar ectoine molecules. 
Furthermore direct stabilizing interactions like intra- and 
intermolecular hydrogen bonds are nearly absent as our 
results have shown. By studying the properties of a 
water-hydroxyectoine mixture we have indicated a sig- 
nificant increase of the excess free energy in combina- 
tion with decreased excess entropy values. The observed 
entropy deficit can be validated by an increase of re- 
laxational orientation times, increase of water hydrogen 
bond lifetimes as well as a decreased diffusivity. The free 
energy increase can be also related to a lower average 
number of hydrogen bonds between water molecules. 
By a combination of these observations, we propose the 
following mechanism which is responsible for the preser- 
vation of a proteins compact configuration. Direct inter- 
actions between hydroxyectoine and the macromolecule 
can be neglected due to the absence of intermolecular 
hydrogen bonds. Due to the increased values for the po- 
lar molecules free solvation energy, it can be assumed 
that a swelling of hydrophilic surface areas is less ener- 
getically favorable in presence of a highly concentrated 
solution of hydroxyectoine compared to a pure water so- 
lution. This becomes also obvious by the values for the 
hydrophilic/hydrophobic solvent accessible surface area 
for CI2 shown in Tab. HI It is remarkable that the size of 
the hydrophobic area increases in presence of hydroxyec- 
toine whereas the hydrophilic area decreases. 
As a direct consequence, apolar hydrophobic molecules 
have decreased energy values for large hydroxyectoine 
concentrations as it has been validated for IIO-CII2-OII 
(data not shown). 

This effect is mainly responsible for the stabilization of 
the proteins compact state. It is widely accepted that 
protein folding to the compact state is initiated by a hy- 



drophobic collapse [66|, |68| where hydrophobic residues 
are shielded in the interior whereas polar residues stay in 
direct contact with the solvent [g^ [g^ . 
Considering the reverse direction, unfolding is correlated 
with a swelling of the total solvent accessible surface area. 
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In the presence of hydroxyectoine, an increase of the hy- 
drophihc surface area is energetically unfavorable. Com- 
pared to a pure solution, the compact native state is sta- 
bilized due to an increase of the free energy barriers such 
that swelling is suppressed. 

It has to be noticed that the proposed mechanism is uni- 
versal for all proteins with well-defined hydrophilic sur- 
face areas regardless of the chemical details. 
The deficits for the free solvation energy in presence of 
compatible solutes are in good agreement to experimen- 
tal results i35i] and to the conclusions of the transfer free 
energy model [2^ which assumes a decrease of solvent 
quality for large co-solvent concentrations. Hence the in- 
teraction of the macromolecules with the solvent becomes 
less favorable such that a shrinkage process due to solvo- 
phobic effects can be observed. Detailed deviations 
from the transfer free energy model in our simulation 
are given by a negligible excess number of intramolecular 
hydrogen bonds to stabilize the proteins compact struc- 
ture. 

Our results validate that solvophobic effects are strong 
enough to stabilize the preservation of the compact native 
state as it can be also seen by the free energy landscapes 
shown in Fig. |4] and as it has been proposed in Ref. [35j . 
Unfolding is energetically less preferable in presence of 
hydroxyectoine such that the equilibrium constant even 
at harsh conditions is shifted towards the native struc- 
ture. These results are in good agreement to experimen- 
tal findings [s^. 

Finally we are also able to explain the results of a recent 
publication fld\, where it has been found that the melt- 
ing temperature for Met-Enkephalin remains unchanged 
in presence of ectoine. By detailed analysis of Met- 
Enkephalin, it comes out that the the ratio of the hy- 
drophobic surface area to the total solvent accessible sur- 
face area is larger compared to CI2. Due to the discussion 
above, we conclude that the stabilizing effects of the co- 
solvent on the hydrophilic regions are less important for 
Met-Enkephalin. In addition it has been also validated 
that a well-defined unfolding process for Met-Enkephalin 
is hard to observe due to low energy barriers between 
different conformations [H^ . 

Summarizing all results, we propose a simple mechanism 
in close agreement to the transfer free energy model to 
explain the stabilizing effect of extremolytes on the native 
state of proteins. Our results imply that the solvophobic 
effect is the main driving force to shift the equilibrium 
constant to the native form in absence of direct intra- and 
intcrmolccular interactions. The strategy of extremolyte 
production in microbacteriae under environmental stress 
can be explained by the proposed mechanism which is 
validated by significant variations of the solvent proper- 
ties. 
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